Efficient design of synthetic gene circuits under cell-to-cell variability.

BACKGROUND: Synthetic biologists use and combine diverse biological parts to build systems such as genetic circuits that perform desirable functions in, for example, biomedical or industrial applications. Computer-aided design methods have been developed to help choose appropriate network structures and biological parts for a given design objective. However, they almost always model the behavior of the network in an average cell, despite pervasive cell-to-cell variability. RESULTS: Here, we present a computational framework and an efficient algorithm to guide the design of synthetic biological circuits while accounting for cell-to-cell variability explicitly. Our design method integrates a Non-linear Mixed-Effects (NLME) framework into a Markov Chain Monte-Carlo (MCMC) algorithm for design based on ordinary differential equation (ODE) models. The analysis of a recently developed transcriptional controller demonstrates first insights into design guidelines when trying to achieve reliable performance under cell-to-cell variability. CONCLUSION: We anticipate that our method not only facilitates the rational design of synthetic networks under cell-to-cell variability, but also enables novel applications by supporting design objectives that specify the desired behavior of cell populations.

A quantitative model for virus uncoating predicts influenza A infectivity.

For virus infection of new host cells, the disassembly of the protective outer protein shell (capsid) is a critical step, but the mechanisms and host-virus interactions underlying the dynamic, active, and regulated uncoating process are largely unknown. Here, we develop an experimentally supported, multiscale kinetics model that elucidates mechanisms of influenza A virus (IAV) uncoating in cells. Biophysical modeling demonstrates that interactions between capsid M1 proteins, host histone deacetylase 6 (HDAC6), and molecular motors can physically break the capsid in a tug-of-war mechanism. Biochemical analysis and biochemical-biophysical modeling identify unanchored ubiquitin chains as essential and allow robust prediction of uncoating efficiency in cells. Remarkably, the different infectivity of two clinical strains can be ascribed to a single amino acid variation in M1 that affects binding to HDAC6. By identifying crucial modules of viral infection kinetics, the mechanisms and models presented here could help formulate novel strategies for broad-range antiviral treatment.

Assumptions on decision making and environment can yield multiple steady states in microbial community models.

BACKGROUND: Microbial community simulations using genome scale metabolic networks (GSMs) are relevant for many application areas, such as the analysis of the human microbiome. Such simulations rely on assumptions about the culturing environment, affecting if the culture may reach a metabolically stationary state with constant microbial concentrations. They also require assumptions on decision making by the microbes: metabolic strategies can be in the interest of individual community members or of the whole community. However, the impact of such common assumptions on community simulation results has not been investigated systematically. RESULTS: Here, we investigate four combinations of assumptions, elucidate how they are applied in literature, provide novel mathematical formulations for their simulation, and show how the resulting predictions differ qualitatively. Our results stress that different assumption combinations give qualitatively different predictions on microbial coexistence by differential substrate utilization. This fundamental mechanism is critically under explored in the steady state GSM literature with its strong focus on coexistence states due to crossfeeding (division of labor). Furthermore, investigating a realistic synthetic community, where the two involved strains exhibit no growth in isolation, but grow as a community, we predict multiple modes of cooperation, even without an explicit cooperation mechanism. CONCLUSIONS: Steady state GSM modelling of microbial communities relies both on assumed decision making principles and environmental assumptions. In principle, dynamic flux balance analysis addresses both. In practice, our methods that address the steady state directly may be preferable, especially if the community is expected to display multiple steady states.

The motor domain of the kinesin Kip2 promotes microtubule polymerization at microtubule tips.

Kinesins are microtubule-dependent motor proteins, some of which moonlight as microtubule polymerases, such as the yeast protein Kip2. Here, we show that the CLIP-170 ortholog Bik1 stabilizes Kip2 at microtubule ends where the motor domain of Kip2 promotes microtubule polymerization. Live-cell imaging and mathematical estimation of Kip2 dynamics reveal that disrupting the Kip2-Bik1 interaction aborts Kip2 dwelling at microtubule ends and abrogates its microtubule polymerization activity. Structural modeling and biochemical experiments identify a patch of positively charged residues that enables the motor domain to bind free tubulin dimers alternatively to the microtubule shaft. Neutralizing this patch abolished the ability of Kip2 to promote microtubule growth both in vivo and in vitro without affecting its ability to walk along microtubules. Our studies suggest that Kip2 utilizes Bik1 as a cofactor to track microtubule tips, where its motor domain then recruits free tubulin and catalyzes microtubule assembly.

Functional comparison of metabolic networks across species.

Metabolic phenotypes are pivotal for many areas, but disentangling how evolutionary history and environmental adaptation shape these phenotypes is an open problem. Especially for microbes, which are metabolically diverse and often interact in complex communities, few phenotypes can be determined directly. Instead, potential phenotypes are commonly inferred from genomic information, and rarely were model-predicted phenotypes employed beyond the species level. Here, we propose sensitivity correlations to quantify similarity of predicted metabolic network responses to perturbations, and thereby link genotype and environment to phenotype. We show that these correlations provide a consistent functional complement to genomic information by capturing how network context shapes gene function. This enables, for example, phylogenetic inference across all domains of life at the organism level. For 245 bacterial species, we identify conserved and variable metabolic functions, elucidate the quantitative impact of evolutionary history and ecological niche on these functions, and generate hypotheses on associated metabolic phenotypes. We expect our framework for the joint interpretation of metabolic phenotypes, evolution, and environment to help guide future empirical studies.

Multivalency ensures persistence of a +TIP body at specialized microtubule ends.

Microtubule plus-end tracking proteins (+TIPs) control microtubule specialization and are as such essential for cell division and morphogenesis. Here we investigated interactions and functions of the budding yeast Kar9 network consisting of the core +TIP proteins Kar9 (functional homologue of APC, MACF and SLAIN), Bim1 (orthologous to EB1) and Bik1 (orthologous to CLIP-170). A multivalent web of redundant interactions links the three +TIPs together to form a '+TIP body' at the end of chosen microtubules. This body behaves as a liquid condensate that allows it to persist on both growing and shrinking microtubule ends, and to function as a mechanical coupling device between microtubules and actin cables. Our study identifies nanometre-scale condensates as effective cellular structures and underlines the power of dissecting the web of low-affinity interactions driving liquid-liquid phase separation in order to establish how condensation processes support cell function.

Cellular and computational models reveal environmental and metabolic interactions in MMUT-type methylmalonic aciduria.

Methylmalonyl-coenzyme A (CoA) mutase (MMUT)-type methylmalonic aciduria is a rare inherited metabolic disease caused by the loss of function of the MMUT enzyme. Patients develop symptoms resembling those of primary mitochondrial disorders, but the underlying causes of mitochondrial dysfunction remain unclear. Here, we examined environmental and genetic interactions in MMUT deficiency using a combination of computational modeling and cellular models to decipher pathways interacting with MMUT. Immortalized fibroblast (hTERT BJ5ta) MMUT-KO (MUTKO) clones displayed a mild mitochondrial impairment in standard glucose-based medium, but they did not to show increased reliance on respiratory metabolism nor reduced growth or viability. Consistently, our modeling predicted MUTKO specific growth phenotypes only for lower extracellular glutamine concentrations. Indeed, two of three MMUT-deficient BJ5ta cell lines showed a reduced viability in glutamine-free medium. Further, growth on 183 different carbon and nitrogen substrates identified increased NADH (nicotinamide adenine dinucleotide) metabolism of BJ5ta and HEK293 MUTKO cells compared with controls on purine- and glutamine-based substrates. With this knowledge, our modeling predicted 13 reactions interacting with MMUT that potentiate an effect on growth, primarily those of secondary oxidation of propionyl-CoA, oxidative phosphorylation and oxygen diffusion. Of these, we validated 3-hydroxyisobutytyl-CoA hydrolase (HIBCH) in the secondary propionyl-CoA oxidation pathway. Altogether, these results suggest compensation for the loss of MMUT function by increasing anaplerosis through glutamine or by diverting flux away from MMUT through the secondary propionyl-CoA oxidation pathway, which may have therapeutic relevance.

Cell region fingerprints enable highly precise single-cell tracking and lineage reconstruction.

Experimental studies of cell growth, inheritance and their associated processes by microscopy require accurate single-cell observations of sufficient duration to reconstruct the genealogy. However, cell tracking-assigning identical cells on consecutive images to a track-is often challenging, resulting in laborious manual verification. Here, we propose fingerprints to identify problematic assignments rapidly. A fingerprint distance compares the structural information contained in the low frequencies of a Fourier transform to measure the similarity between cells in two consecutive images. We show that fingerprints are broadly applicable across cell types and image modalities, provided the image has sufficient structural information. Our tracker (TracX) uses fingerprints to reject unlikely assignments, thereby increasing tracking performance on published and newly generated long-term data sets. For Saccharomyces cerevisiae, we propose a comprehensive model for cell size control at the single-cell and population level centered on the Whi5 regulator, demonstrating how precise tracking can help uncover previously undescribed single-cell biology.

Model-based inference of neutralizing antibody avidities against influenza virus.

To assess the response to vaccination, quantity (concentration) and quality (avidity) of neutralizing antibodies are the most important parameters. Specifically, an increase in avidity indicates germinal center formation, which is required for establishing long-term protection. For influenza, the classical hemagglutination inhibition (HI) assay, however, quantifies a combination of both, and to separately determine avidity requires high experimental effort. We developed from first principles a biophysical model of hemagglutination inhibition to infer IgG antibody avidities from measured HI titers and IgG concentrations. The model accurately describes the relationship between neutralizing antibody concentration/avidity and HI titer, and explains quantitative aspects of the HI assay, such as robustness to pipetting errors and detection limit. We applied our model to infer avidities against the pandemic 2009 H1N1 influenza virus in vaccinated patients (n = 45) after hematopoietic stem cell transplantation (HSCT) and validated our results with independent avidity measurements using an enzyme-linked immunosorbent assay with urea elution. Avidities inferred by the model correlated with experimentally determined avidities (ρ = 0.54, 95% CI = [0.31, 0.70], P < 10-4). The model predicted that increases in IgG concentration mainly contribute to the observed HI titer increases in HSCT patients and that immunosuppressive treatment is associated with lower baseline avidities. Since our approach requires only easy-to-establish measurements as input, we anticipate that it will help to disentangle causes for poor vaccination outcomes also in larger patient populations. This study demonstrates that biophysical modelling can provide quantitative insights into agglutination assays and complement experimental measurements to refine antibody response analyses.

PolyRound: polytope rounding for random sampling in metabolic networks.

SUMMARY: Random flux sampling is a powerful tool for the constraint-based analysis of metabolic networks. The most efficient sampling method relies on a rounding transform of the constraint polytope, but no available rounding implementation can round all relevant models. By removing redundant polytope constraints on the go, PolyRound simplifies the numerical problem and rounds all the 108 models in the BiGG database without parameter tuning, compared to ∼50% for the state-of-the-art implementation. AVAILABILITY AND IMPLEMENTATION: The implementation is available on gitlab: https://gitlab.com/csb.ethz/PolyRound. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Association of Host Factors With Antibody Response to Seasonal Influenza Vaccination in Allogeneic Hematopoietic Stem Cell Transplant Patients.

BACKGROUND: Influenza vaccination efficacy is reduced after hematopoietic stem cell transplantation (HSCT) and patient factors determining vaccination outcomes are still poorly understood. METHODS: We investigated the antibody response to seasonal influenza vaccination in 135 HSCT patients and 69 healthy volunteers (HVs) in a prospective observational multicenter cohort study. We identified patient factors associated with hemagglutination inhibition titers against A/California/2009/H1N1, A/Texas/2012/H3N2, and B/Massachusetts/2012 by multivariable regression on the observed titer levels and on seroconversion/seroprotection categories for comparison. RESULTS: Both regression approaches yielded consistent results but regression on titers estimated associations with higher precision. HSCT patients required 2 vaccine doses to achieve average responses comparable to a single dose in HVs. Prevaccination titers were positively associated with time after transplantation, confirming that HSCT patients can elicit potent antibody responses. However, an unrelated donor, absolute lymphocyte counts below the normal range, and treatment with calcineurin inhibitors lowered the odds of responding. CONCLUSIONS: HSCT patients show a highly heterogeneous vaccine response but, overall, patients benefited from the booster shot and can acquire seroprotective antibodies over the years after transplantation. Several common patient factors lower the odds of responding, urging identification of additional preventive strategies in the poorly responding groups. CLINICAL TRIALS REGISTRATION: NCT03467074.

Probabilistic thermodynamic analysis of metabolic networks.

MOTIVATION: Random sampling of metabolic fluxes can provide a comprehensive description of the capabilities of a metabolic network. However, current sampling approaches do not model thermodynamics explicitly, leading to inaccurate predictions of an organism's potential or actual metabolic operations. RESULTS: We present a probabilistic framework combining thermodynamic quantities with steady-state flux constraints to analyze the properties of a metabolic network. It includes methods for probabilistic metabolic optimization and for joint sampling of thermodynamic and flux spaces. Applied to a model of Escherichia coli, we use the methods to reveal known and novel mechanisms of substrate channeling, and to accurately predict reaction directions and metabolite concentrations. Interestingly, predicted flux distributions are multimodal, leading to discrete hypotheses on E.coli's metabolic capabilities. AVAILABILITY AND IMPLEMENTATION: Python and MATLAB packages available at https://gitlab.com/csb.ethz/pta. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A rationally engineered decoder of transient intracellular signals.

Cells can encode information about their environment by modulating signaling dynamics and responding accordingly. Yet, the mechanisms cells use to decode these dynamics remain unknown when cells respond exclusively to transient signals. Here, we approach design principles underlying such decoding by rationally engineering a synthetic short-pulse decoder in budding yeast. A computational method for rapid prototyping, TopoDesign, allowed us to explore 4122 possible circuit architectures, design targeted experiments, and then rationally select a single circuit for implementation. This circuit demonstrates short-pulse decoding through incoherent feedforward and positive feedback. We predict incoherent feedforward to be essential for decoding transient signals, thereby complementing proposed design principles of temporal filtering, the ability to respond to sustained signals, but not to transient signals. More generally, we anticipate TopoDesign to help designing other synthetic circuits with non-intuitive dynamics, simply by assembling available biological components.

Model-based integration of genomics and metabolomics reveals SNP functionality in Mycobacterium tuberculosis.

Human tuberculosis is caused by members of the Mycobacterium tuberculosis complex (MTBC) that vary in virulence and transmissibility. While genome-wide association studies have uncovered several mutations conferring drug resistance, much less is known about the factors underlying other bacterial phenotypes. Variation in the outcome of tuberculosis infection and diseases has been attributed primarily to patient and environmental factors, but recent evidence indicates an additional role for the genetic diversity among MTBC clinical strains. Here, we used metabolomics to unravel the effect of genetic variation on the strain-specific metabolic adaptive capacity and vulnerability. To define the functionality of single-nucleotide polymorphisms (SNPs) systematically, we developed a constraint-based approach that integrates metabolomic and genomic data. Our model-based predictions correctly classify SNP effects in pyruvate kinase and suggest a genetic basis for strain-specific inherent baseline susceptibility to the antibiotic para-aminosalicylic acid. Our method is broadly applicable across microbial life, opening possibilities for the development of more selective treatment strategies.

Efficient manipulation and generation of Kirchhoff polynomials for the analysis of non-equilibrium biochemical reaction networks.

Kirchhoff polynomials are central for deriving symbolic steady-state expressions of models whose dynamics are governed by linear diffusion on graphs. In biology, such models have been unified under a common linear framework subsuming studies across areas such as enzyme kinetics, G-protein coupled receptors, ion channels and gene regulation. Due to 'history dependence' away from thermodynamic equilibrium, these models suffer from a (super) exponential growth in the size of their symbolic steady-state expressions and, respectively, Kirchhoff polynomials. This algebraic explosion has limited applications of the linear framework. However, recent results on the graph-based prime factorization of Kirchhoff polynomials may help subdue the combinatorial complexity. By prime decomposing the graphs contained in an expression of Kirchhoff polynomials and identifying the graphs giving rise to equal polynomials, we formulate a coarse-grained variant of the expression suitable for symbolic simplification. We devise criteria to efficiently test the equality of Kirchhoff polynomials and propose two heuristic algorithms to explicitly generate individual Kirchhoff polynomials in a compressed form; they are inspired by algebraic simplifications but operate on the corresponding graphs. We illustrate the practicality of the developed theory and algorithms for a diverse set of graphs of different sizes and for non-equilibrium gene regulation analyses.

TopoFilter: a MATLAB package for mechanistic model identification in systems biology.

BACKGROUND: To develop mechanistic dynamic models in systems biology, one often needs to identify all (or minimal) representations of the biological processes that are consistent with experimental data, out of a potentially large set of hypothetical mechanisms. However, a simple enumeration of all alternatives becomes quickly intractable when the number of model parameters grows. Selecting appropriate dynamic models out of a large ensemble of models, taking the uncertainty in our biological knowledge and in the experimental data into account, is therefore a key current problem in systems biology. RESULTS: The TopoFilter package addresses this problem in a heuristic and automated fashion by implementing the previously described topological filtering method for Bayesian model selection. It includes a core heuristic for searching the space of submodels of a parametrized model, coupled with a sampling-based exploration of the parameter space. Recent developments of the method allow to balance exhaustiveness and speed of the model space search, to efficiently re-sample parameters, to parallelize the search, and to use custom scoring functions. We use a theoretical example to motivate these features and then demonstrate TopoFilter's applicability for a yeast signaling network with more than 250'000 possible model structures. CONCLUSIONS: TopoFilter is a flexible software framework that makes Bayesian model selection and reduction efficient and scalable to network models of a complexity that represents contemporary problems in, for example, cell signaling. TopoFilter is open-source, available under the GPL-3.0 license at https://gitlab.com/csb.ethz/TopoFilter. It includes installation instructions, a quickstart guide, a description of all package options, and multiple examples.

Controlling cell-to-cell variability with synthetic gene circuits.

Cell-to-cell variability originating, for example, from the intrinsic stochasticity of gene expression, presents challenges for designing synthetic gene circuits that perform robustly. Conversely, synthetic biology approaches are instrumental in uncovering mechanisms underlying variability in natural systems. With a focus on reducing noise in individual genes, the field has established a broad synthetic toolset. This includes noise control by engineering of transcription and translation mechanisms either individually, or in combination to achieve independent regulation of mean expression and its variability. Synthetic feedback circuits use these components to establish more robust operation in closed-loop, either by drawing on, but also by extending traditional engineering concepts. In this perspective, we argue that major conceptual advances will require new theory of control adapted to biology, extensions from single genes to networks, more systematic considerations of origins of variability other than intrinsic noise, and an exploration of how noise shaping, instead of noise reduction, could establish new synthetic functions or help understanding natural functions.

Remote control of microtubule plus-end dynamics and function from the minus-end.

In eukaryotes, the organization and function of the microtubule cytoskeleton depend on the allocation of different roles to individual microtubules. For example, many asymmetrically dividing cells differentially specify microtubule behavior at old and new centrosomes. Here we show that yeast spindle pole bodies (SPBs, yeast centrosomes) differentially control the plus-end dynamics and cargoes of their astral microtubules, remotely from the minus-end. The old SPB recruits the kinesin motor protein Kip2, which then translocates to the plus-end of the emanating microtubules, promotes their extension and delivers dynein into the bud. Kip2 recruitment at the SPB depends on Bub2 and Bfa1, and phosphorylation of cytoplasmic Kip2 prevents random lattice binding. Releasing Kip2 of its control by SPBs equalizes its distribution, the length of microtubules and dynein distribution between the mother cell and its bud. These observations reveal that microtubule organizing centers use minus to plus-end directed remote control to individualize microtubule function.

Publisher Correction: Microbial network disturbances in relapsing refractory Crohn's disease.

Owing to an error during typesetting, a number of references were deleted from the Methods reference list. This altered all of the references in the Methods section and some of the references in Extended Data Fig. 5, making them inaccurate. References 121-134 were added back into the Methods reference list, and the references in the Methods section and in Extended Data Fig. 5 were renumbered accordingly. The error has been corrected in the PDF and HTML versions of this article.

Microbial network disturbances in relapsing refractory Crohn's disease.

Inflammatory bowel diseases (IBD) can be broadly divided into Crohn's disease (CD) and ulcerative colitis (UC) from their clinical phenotypes. Over 150 host susceptibility genes have been described, although most overlap between CD, UC and their subtypes, and they do not adequately account for the overall incidence or the highly variable severity of disease. Replicating key findings between two long-term IBD cohorts, we have defined distinct networks of taxa associations within intestinal biopsies of CD and UC patients. Disturbances in an association network containing taxa of the Lachnospiraceae and Ruminococcaceae families, typically producing short chain fatty acids, characterize frequently relapsing disease and poor responses to treatment with anti-TNF-α therapeutic antibodies. Alterations of taxa within this network also characterize risk of later disease recurrence of patients in remission after the active inflamed segment of CD has been surgically removed.